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Abstract. There is only one fully supported ergodic invariant probability 
measure for the adic transformation on the space of infinite paths in the graph 
that underlies the Eulerian numbers. This result may partially justify a fre- 
quent assumption about the equidistribution of random permutations. 



1. Introduction 

Wc give a new proof of the main result of [?] , crgodicity of the symmetric measure 
for the Euler adic transformation. In fact we prove a stronger result by a differ- 
ent method: the symmetric measure is the only fully supported ergodic invariant 
probability measure for the Euler adic system. 

This result has a probabilistic interpretation. For each n = 0, 1,2, ...,k — 
0, 1, ... n, the Eulerian number A(n, k) is the number of permutations i\ii . . . z ra +i 
of {1, . . . , n + 1} with exactly k rises (indices j = 1,2, ... ,n with ij < ij+i) and 
n — k falls (indices j = 1,2, ...,n with ij > ij+i). (See, for example, [?] for 
basic facts about the Eulerian numbers.) Besides their combinatorial importance, 
these numbers are also of interest in connection with the statistics of rankings: see, 
for example, [?],[?],[?],[?], and [?]. In studying random permutations, it is often 
assumed that all permutations are equally likely, each permutation of length n + 1 
occurring with probability l/(n + 1)!. Our main result implies that in a sense 
(see Remark[T]) there is a unique way to choose permutations at random so that (1) 
consistency holds: the distribution on permutations of {1, ... , n + m} induces, upon 
deletion of {n + 1, . . . , to}, the distribution on {1, . . . , n}; (2) any two permutations 
of the same length which have the same number of rises are equally likely; and 
(3) every permutation has positive probability. This unique way is to make all 
permutations of the same length, no matter the number of rises, equally likely. 

This conclusion is achieved by means of ergodic theory, studying invariant mea- 
sures for the adic (Bratteli-Vershik) transformation on the space of infinite paths 
on the Euler graph. We do not see how to prove this result by a straightforward 
combinatorial argument, because it is very difficult to make comparative asymp- 
totic estimates of Eulerian numbers. Indeed, the key to our argument is to finesse 
this difficulty (in the proof of Proposition [3]) by setting up a situation in which the 
dominant terms of the two expressions being compared are identical. 

Not having the space here to survey all the literature on adic systems, we remark 
just that in certain senses they represent all dynamical systems, also have strong 
connections with group representations, C* algebras, probability theory, and com- 
binatorics, and thus facilitate the study of important problems in those areas. See 
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[?, ?, ?, ?, ?, ?] and the references that these works contain for background. Basic 
attributes of adic systems are the various kinds of dimension groups, which provide 
invariants, and the states or traces, which correspond to adic-invariant, or "central" , 
measures — see, for example, [?]. These have been determined for many stationary 
and essentially simple adic systems and some classes of nonstationary and nonsim- 
ple ones [?, ?], but not yet for many key individual examples. For the Pascal adic, 
the identification of the ergodic invariant measures as the one-parameter family of 
Bernoulli measures on the 2-shift is closely connected with de Finetti's Theorem, 
the Hewitt-Savage 0,1 Law, and uniform distribution of the Kakutani interval split- 
ting procedure — see the discussion in [?] . Identifying the ergodic measures for new 
natural examples such as the Euler adic will also have consequences in other ar- 
eas: mentioned below are connections with random permutations (Remark [T]) and 
reinforced random walks (Remark [4]). 

2. The Euler adic system 

The Euler graph T is a Bratteli diagram with levels n = 0, 1, 2, . . . . There are 
?i + 1 vertices of each level n, labelled (n,k), < k < n, and each vertex (n, k) is 
connected to vertex {n + 1, k) by k + 1 edges and to vertex (n + 1, k + 1) by n — k + 1 
edges. The number of paths into any vertex (n, k) is the Eulerian number A(n, k). 
The Eulerian numbers have the recursion 

(1) A(n, k) = (n-k + l)A(n - 1, k - 1) + (k + l)A(n - 1, k), 

where A(0, 0) = 1 and by convention A(n, k) = for k ^ {0, 1, n}. 

Level 




Figure 1 . The Euler graph down to level 5 



Define X to be the space of infinite edge paths in this graph: X = {x — (x n )\n = 
0, 1, 2, . . . }, each x n being an edge from level n to level n + 1. The vertex through 
which x passes at level n is denoted by (n, k n {x)). We say that an edge is a left turn 
if it connects vertices (n, k) and (n + 1, k) and a right turn if it connects vertices 
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(n, k) and (n + 1, k + 1). X is a metric space in the usual way: d(x, y) = 2~ J , where 
j = inf{i|x; 7^ yi}. The cylinder sets, where a finite number of edges are specified, 
are clopen sets that generate the topology of X. 

As with other Bratteli diagrams, we define a partial order on the edges in the 
diagram which extends to the entire space of paths. Two edges are comparable if 
they terminate in the same vertex. For each vertex, totally order the set of edges 
that terminate at that vertex. This is pictured by agreeing that the minimal edge 
between two vertices is the leftmost edge, and the edges increase from left to right. 
Two paths, x, y, are comparable if they agree after some level. Then x is less than 
y if the last edges of x and y that do not agree, Xi and yi, are such that Xi is less 
than yi in the edge ordering. We define the adic transformation T on the set of 
non-maximal paths in X to map a path to the next largest path according to this 
induced partial order. Then two paths arc in the same orbit if and only if they are 
comparable. Figure [2] shows a particular path x and its image y — Tx. There are 
a countable number of paths for which there is no next largest path, we denote the 
set of such paths by X max . In particular for each k in N U {oo}, we define £ ma x(fc) 
to be the infinite edge path that follows the unique finite path from (0,0) to (k, k) 
and for n > k, follows the right most edge connecting vertices (n, k) and (n + 1, k). 
In a similar fashion there are a countable number of paths for which there is no 
next smaller path, denoted by X m i n , and for every k in N U {oo} a unique path 
^min(fc) S ^min- The transformation T maps x maxfe to x min(k) . 




Figure 2. The dotted path maps to the dashed path 

The symmetric measure, rj, on the infinite path space X of the Euler graph is the 
Borel probability measure that for each n gives every cylinder of length n starting 
at the root vertex the same measure. Clearly rj is T-invariant. The measure of any 
cylinder set can be computed by multiplying weights on the edges, each weight on 
an edge connecting level n to level n+l being l/(n + 2). We can think of the weights 
as assigning equal probabilities to all the allowed steps for a random walker who 
starts at the root and descends step by step to form an infinite path x G X. The 
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Figure 3. The symmetric measure given by weights on edges 



main result of [?], proved by probabilistic methods, is that the symmetric measure 
is ergodic. 

There is a bijective correspondence between paths (or cylinders) of length no 
starting at the root vertex and terminating at vertex (no,ko) and permutations 
of {1, 2, . . . , no + 1} with ko rises. This correspondence is crucial for our proof of 
Proposition [3l which is the essential component of our main result. Consider the 
cylinder set defined by the single edge connecting the vertex (0, 0) to the vertex 
(1,0). This cylinder set is of length 1 with 1 left turn, and we assign to it the 
permutation 21, which has one fall. Likewise, the cylinder set defined by the single 
edge connecting the vertex (0,0) to the vertex (1, 1) is of length 1 with one right 
turn, and we assign to it the permutation 12, which has one rise. When a cylinder F 
of length n, corresponding to the permutation ir(F) of {1, 2, . . . , n + 1}, is extended 
by an edge from level n to level n + 1, we extend ir(F) in a unique way to a 
permutation of {1, 2, . . . ,n + 2}, as follows. If F is extended by a left turn down 
the i'th edge connecting (n, k) to (n + 1, k), insert n + 2 into tt(F) in the i'th place 
that adds an additional fall to the total number of falls of tt(F). Likewise, if F 
is extended by a right turn down the i'th edge connecting (n, k) to {n + 1, k + 1), 
insert n + 2 into ir(F) in the i'th place that adds an additional rise to the total 
number of rises of tt(F). 

Example 1. Wc shall determine which cylinder F corresponds to the permutation 
tt(F) = 2341. When we delete 34, we have the permutation 21, which means the 
first edge of F is the unique edge connecting vertex (0,0) to (1,0). When the 3 is 
added, it is added into the middle of 21, creating the permutation 231, which has 
one rise and one fall. Therefore a rise has been added, and there is no earlier space 
that the 3 could have been placed to add a rise to 21. Therefore, the second edge 
of F is the first edge connecting vertex (1,0) to vertex (2,1). Since ir(F) has 2 rises 
and one fall, the addition of 4 into 2341 adds another rise. It is added in the first 
place it could be added in order to add a rise, meaning that the third edge of F 
must be the first edge connecting vertex (2,1) to (3,2). See Figure |U 

This correspondence produces a labeling of infinite paths in the Euler graph 
starting at the root; then the path space X corresponds to the set of all linear 
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orderings of N = {1,2,3,...}, and the adic transformation T can be thought of as 
moving from an ordering to its successor. 




3. Proof of the main result — uniqueness of the symmetric measure 

The dimension dim(n, k) of a vertex (n, k) is defined to be the number of paths 
connecting the root vertex (0,0) to the vertex (n, k). For any cylinder F define 
dim(F, (n, k)) to be the number of paths in F that connect the root vertex (0, 0) 
to the vertex (n, k). If F is a cylinder starting at the root vertex (0,0), then 
dim(_F, (n, k)) is the number of paths from the terminal vertex of F to (n, k). Denote 
by I the er-algcbra of T-invariant Borel measurable subsets of X. The following is 
a well-known result about the measures of cylinders in an adic situation. 

Lemma 1 (Vershik [?],[?]). Let fi be an invariant probability measure for the Euler 
adic transformation. Then for every cylinder set F and [i-almost every x G X , 

Hm dim(F,(n,k n (x))) = ^ (y f \X). 
n->oo dim((n, k n (x))) 

It is clear from the Markov property and the patterns of weights on the edges 
that with respect to n almost every path x £ X has infinitely many left and right 
turns. 

Lemma 2. Let fi be an invariant fully- supported ergodic probability measure for 
the Euler adic transformation. For [i-almost every x G X , there are infinitely many 
left and right turns (i.e., k n (x) and n — k n (x) are unbounded a.e.). 

Proof. For each K = 1, 2, . . . , let 

Ak = {% G X\x has no more than K right turns}. 

Then Ak is a proper closed T-invariant set. Since fi is ergodic and fully supported, 
h(Ak ) = 0. Similarly, the measure of the set of paths with a bounded number of 
left turns is 0. □ 

Proposition. // fi is an invariant probability measure for the Euler adic transfor- 
mation such that k n (x) and n — k n {x) are unbounded a.e., then \i = 77. 
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Proof. For any string w on an ordered alphabet denote by r(w) the number of 
rises in w and by f(w) the number of falls in w (defined as above). Let F and 
F' be cylinder sets in X specified by fixing the first no edges, and let tt(F) and 
tt(F') be the permutations assigned to them by the correspondence described in the 
preceding section. Suppose that the paths corresponding to F and F' terminate in 
the vertices (no, fco) an d (no, k ) respectively. Fix n ^> no and fc 3> fco- 

We will show that p(F) = p(F'), and hence p = n. If fco = k' , then p(F) = p(F') 
by the T-invariance of p, so there is nothing to prove. We will deal with the case 
fc ^ k' Q . 

Example 2. Let F be the dotted cylinder, and F' the dashed. Then tt(F) = 213 
and tt(F') = 132. For purposes of following an example through, we will let (n, fc) = 
(8,4). 



« 




(n,fc) = (8,4) 



Each path s in T from (no, fco) to (n, fc) corresponds to a permutation a s of 
{1, 2, . . . , n + 1} with fc rises in which 1, 2, . . . , no + 1 appear in the order tt(F). 
Counting dim(i^, (n, fc)) is equivalent to counting the number of distinct such a s - 
Each such permutation a s has associated to it a permutation t(a s ) of {no+2, . . . , n+ 
1} obtained by deleting 1, 2, ... , no + 1 from cr s , see Example [31 Taking a reverse 
view, one obtains er s from p = t(o~ s ) by inserting 1, 2, . . . , no + 1 from left to right, 
in the order prescribed by tv(F), into p. 

We define a cluster in o~ s to be a subset of {1, 2, . . . , n + 1} whose members 
are found consecutively in a si with no elements of {no + 2, . . . , n + 1} separating 
them, in the order prescribed by tt(F). The set M s of clusters in a s is an ordered 
partition of the permutation tt(F), and we define 

r(M s ) = ]T r(c). 

ceM s 

In general, 1 < \M S \ < no + 1 and < r(M s ) < fco. 
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Example 3. n(F) =213 



a si = 297146385 
t(er Si ) = 974685 
M S1 = {2,1,3} 
r(M Sl ) = 



cr s „ = 962471358 
t(a S2 ) = 964758 
M S2 = {2, 13} 
r(M S2 ) = 1 



Given a permutation p of {no + 2, . . . ,n + 1}, < m < |p| + 1, and an ordered 
partition M of 7r(F) with |M| = m, there are C(|p| + 1, m) (the binomial coefficient) 
choices for how to insert the members of M as clusters into the permutation p in 
order to form a permutation a s . But not all of these choices yield valid permutations 
a s , which have exactly fc rises. Looking more closely, we see that placing a cluster 
c G M at the tail end of p or into a rise in p produces a permutation 7j whose 
number of rises is r(p) = r(p) + r(c), while placing c at the beginning or into a fall 
produces p with r(p) = r(p) + r(c) + 1. So we must have 



In order to count the number of valid ways to place the members of M into p, we 
will first determine how many ways there are to place the clusters that will create 
a new rise. There are n — no — (fc — fco) + 1 possible places in p to place a cluster of 
M and create a new rise (the n — uq — (fc — ko) falls and at the end of p), and since 
{c G M\c is placed into a fall or at the beginning of p}\ is k — r(p) — r(M) we must 
choose l(r(p),AI) = k — r(p) — r(M) of them. For each of these possibilities we 
must then choose places to place the remaining clusters of M. There are k — fco + 1 
places, and m — (k — r(p) — r{M)) remaining clusters. Therefore the number of 
ways (given m, M, and r = r(p)) to place the members of M into p in such a way 
as to form a valid permutation cr s , with fc rises and n — k falls, is 

C(n - n - (fc - fc ) + 1, fc - r(p) - r(M)) C(k - k + 1, m - (fc - r(p) - r(M))). 

For each m = l,2,...,7io + l denote by P m (F) the set of ordered partitions M 
of tt(F) such that \M\ = m. For each r = 0, 1, . . . , n — no — 1 denote by Q(n, no, r) 
the set of permutations of {no + 2, . . . , n + 1} with exactly r rises, so that 

\Q(n,n ,r)\ = A(n - n - l,r). 
Example 4. tt(F) = 213, (n, fc) = (8,4). 



In order to compute dim(i 7 ', (n, fc)) we will partition the set of permutations <j s 
of {1, 2, which have exactly fc rises and a subpermutation ir{F) in the 
following manner. First, partition the set of all such permutations according to 
the cardinality m G {1, . . . , no + 1} of the set of clusters. Now partition further 
by grouping the a s by their corresponding M s G P m (F). Partition yet again by 
grouping the a s by the number of rises, r, in their corresponding permutations t(a s ) 
of {no + 2, . . . , n + 1}. Since the total number of rises in each a s must be fc, and 
there are r{M) rises in M, the minimum such r is fc — r(M) — m and the maximum 



k = r(a s ) =r(p)+r(M)+ 
#{c G M\c is placed into a fall or at the beginning of p}. 



P 1 (F) = {{213}} 
P 2 (F) = {{2,13},{21,3}} 
P 3 (F) = {{2, 1,3}} 
0(8,2,0) = {987654} 
0(8,2,5) = {456789} 



O(8,2,0)| =A(5,0) = 1 
0(8,2,5)| =A{5, 5) = 1 
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is k — r(M). Finally, partition these groups by the distinct permutations t(a s ) of 
{no + 2, . . . , n + 1} which meet the foregoing criteria. Then for each r there are 
\Q(n, Ho,t)| = A(n — no — l,r) of these groups (the same number for each F, to, M), 
and each group has cardinality 

C(n -n -(k-k ) + l,k-r- r(M)) C(k - k + 1, m - (k - r - r(M))). 

Example 5. Fix ir(F) = 213 and (n, k) = (8,4). 

1. Partition by to. 

m = 1 : <t s = . . . 213 . . . 

m = 2 : a s = ...2... 13... or a s = ...21... 3... 
m = 3 : er s = . . . 2 . . . 1 . . . 3 . . . 

2. Fix m = 2, partition by M. 

... 2 ... 13 .. . appears in each <t s 
... 21 ... 3 .. . appears in each a s 

3. Fix m = 2, M = {2, 13}, then r(M) = 1, partition over r{t{a s )). 

r(t(a s )) = 1 

r(t(<r,)) = 2 
r(t(<r a )) = 3 

4. Fix m = 2, M = {2, 13}, r(M) = 1, r(t(a s )) = 3, then partition over t(a s ) 
in Q(8, 2,3). There are A(5,3) = 302 sets in this partition, and each one contains 
C(4,0)C(4, 2) = 6 elements. (To each permutation of 4,5,6,7,8 we can now add 
only one rise, and r(M) is already 1, so all the clusters must be inserted into rises 
or at the end.) 

For ease of notation, for each r, let I = l(r,M) = inf{|M|,fc — r — r(M)} and 
p = p(M) = k- r(M) - \M\. Then dim(F, (n, k)) = 

n + l fe-r(M) 

^2 A(n-n -l 7 r)C(n-n -(k-k ) + l,l)C(k-k + l 7 m-l). 

m=l MEP m (F) r=p 

Regrouping the terms, rewrite the sum as 

k /no+l \ 

(2) dim(F, (n, k)) = ^ I a(F, r, m) I A{n - n - 1, r), 

r=fc-(Ti + l) \m=l / 

where we have 

a(F,r,m)= V" C(tt, - n - (fc - fc ) + 1, /)C(fc - fc + 1, m - Z). 
a/sp,„(F) 

In each term in this sum, the numerator of the first binomial coefficient factor 
simplifies to having / factors chosen from [n — no — (k — fc ) — m, n — ti — (fc — fco)], 
and the numerator of the second has m — I factors chosen from [k — fco — to, A;] . 
Combining the two, the numerator of the binomial factors in each individual term 
in this sum consists of m factors chosen from [n — k — m, n — k] U [k — m, k] , and the 
product of the denominators is bounded above by (to!) 2 . Since to < uq + 1, each 
of the to factors is comparable to (i.e., between two constant multiples of) either k 
or n — fc. Then we see that 



RANDOM PERMUTATIONS AND ERGODICITY FOR THE EULER ADIC 



np + 1 

f3(F,r) = a(F,r,m) 

m— 1 

is a polynomial in k and n — k of degree no + 1. 

Clearly the dominant term of each /3(F, r) as /c and n — k grow large is the one 
of maximum degree, m = no + 1. This is exactly the term a(i 7 ', r, no + 1), which 
corresponds to the partition M of tt(F) into singletons. Then r{M) = 0, and all the 
elements of {1, . . . , no + 1} arc placed, as singleton clusters, into a permutation p of 
{no + 2, . . . , n + 1} which has exactly r rises. Since every element of {1, . . . , no + 1} 
is less than every element of the set {no + 2, . . . , n}, the order of the permutation 
7r(.F) has no effect on a(F, r, no + 1). The main point is that this term is the same 
for ir(F) as for any other permutation ir(F') of {1, 2, . . . , no + 1}: 

a(F, r, n + 1) = a(F', r, n + 1) 

for all r and all F'. When k — r of {1, . . . , no + 1} are put into falls in p or at the 
beginning, and the rest are put into rises or at the end, no matter which elements 
are placed in which slots we always produce a permutation a s with exactly k rises. 

Let us now consider the ratio dim(.F, (n, k))/ dim(i r " , (n, k)) when n and k are 
both very large. Divide top and bottom by the sum on r of the dominant terms 
(taking maximum degree coefficients in k and n — k for each r), 

k 

cn(F,r, n + l)A(n - n - l,r), 

r=k — (n + l) 

which is the same for F and F' . This shows that the ratio is very near 1 when k 
and n — k are both very large. 

Thus if k n (x), n — k n (x) — ► oo a.e. dp, we have for any two cylinders F and F' 
starting at the root vertex and of the same length that 

E^{X F \T){x) = E^{X F .\T){x) a.e. dp. 

Integrating gives p(F) = p(F'), so that p = ?/. □ 

Theorem. The symmetric measure n is ergodic and is the only fully supported 
invariant ergodic Borel probability measure for the Euler adic transformation. 

Proof. If we show that there is an ergodic measure p which has k n (x) and n — k n (x) 
unbounded a.e., then it will follow from the Proposition that p — n, and hence n is 
ergodic and is the only fully-supported T-invariant ergodic measure on X. 

If an ergodic measure has, say, k n {x) bounded on a set of positive measure, then 
k n {x) is bounded a.e., since each set {x\k n (x) < K} is T-invariant. Let So = 0, and 
for each K = 1,2,... let £k be the set of ergodic measures supported on cither 
{x G X\k n {x) < K for all n} or {x G X\n — k n (x) < K for all n}. If no ergodic 
measure has k n (x) and n — k n (x) unbounded a.e., then the set of ergodic measures 
is 

£ = {J£k- 

K 

Form the ergodic decomposition of iy. 

r. OO n 

77= / edP„(e) =Y. edP n (e). 
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If S is the set of paths x 6 X for which both k n (x) and n — k n (x) are unbounded, 
then, from the remark before Lemma [2 r){S) = 1; but, for each K, e(S) = for 
all e in £k . Hence there is an ergodic measure for which k n (x) and n — k n (x) are 



Remark 1. The connection of this Theorem with the statements made in the 
Abstract and Introduction about random permutations can be seen as follows. As 
noted above, the space X of infinite paths in the Euler graph is in correspondence 
with the set £ of linear orderings of N. A cylinder set in X determined by fixing 
an initial path of length n corresponds, as explained above, to a permutation 7r„+i 
in the group S n+ i of permutations of {1, 2, . . . , n + 1}, and thus to the set C(tt„ + i) 
of all elements of C for which 1, 2, . . . , n + 1 appear in the order specified by 7r n+1 . 
This family of clopen cylinder sets defines a compact metrizablc topology and a 
Borel structure on C. One way to speak about "random permutations" would be 
to give a Borel probability measure on C. A Borel probability measure on X is 
T-invariant if and only if the corresponding measure on C assigns, for each n and 

< k < n, equal measure to all cylinders determined by permutations Ttn+i G S n +i 
which have the same fixed number k of rises. According to the Theorem, the only 
such fully supported measure is the one determined by the symmetric measure rj 
on X, which assigns to each basic set C(iT n+ i), for all 7r rl+ i € S n +i, no matter the 
number of rises in 7r„+i, the same measure, l/(n + 1)!. 

Remark 2. It is easy to find measures for the Euler adic which are not fully 
supported, just by restricting to closed T-invariant sets. For example, if we restrict 
to the subgraph consisting of paths x with vertices (n,k n (x)) with k n (x) = or 

1 for all n, then there are invariant measures determined by the edge weights in 
Figure^ provided that a n +\ = a n /(2 — 2na n ), for example a n = l/2(n + 1). All 
such systems are of finite rank (see [?] for the definition) . Since measures which are 
not fully supported do not involve the full richness of the Euler graph, we regard 
them as less interesting than fully supported ergodic measures. 



unbounded a.e.. 



□ 



4. Concluding remarks 



1 - 4a 4 




Figure 5. The invariant measures when T is restricted to the 
paths x with vertices (n, k n (x)) with k n (x) = or 1 for all n. 
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Remark 3. For dynamic properties of the Euler adic with its symmetric measure 
beyond ergodicity, so far we know that this system is totally ergodic (there are 
no roots of unity among its eigenvalues) [?], has entropy 0, and is loosely Bernoulli 
[?] (sec [?] for the definition). Determining whether the system is weakly or even 
strongly mixing and whether it has infinite rank are important, challenging, open 
questions. 

Level 

1 
2 
3 
4 
5 

k = k = l k = 2 fc = 3 k = 4 k = 5 

Figure 6. The Reverse Euler graph 

Remark 4. The Euler adic system models a negatively reinforced random walk 
on the finite graph consisting of two loops based at a single vertex. In other words, 
each time an individual loop is followed the probability that it will be followed 
at the next stage decreases. The adic system based on the Bratteli diagram in 
Figure El which we call the reverse Euler system, models the positively reinforced 
random walk (each time an individual loop is followed the probability that it will 
be followed at the next stage increases) on this graph. Identification of the ergodic 
measures for such systems and determination of their dynamical properties has 
implications for the properties of the random walks. While the Euler adic has only 
one fully supported ergodic measure as we have shown, there is an infinite (one- 
parameter) family of fully supported ergodic measures for the reverse Euler adic. 
The connection between adic systems of this kind and reinforced random walks will 
be explored in a future publication [?] . 
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